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Abstract 

Through the present paper, the code gamgamZZ is presented, which may be used 
to calculate all possible observables related to the process 77 — > ZZ ', in either the Stan- 
dard Model (SM), or the minimal sypersymmetric standard model (MSSM) with real 
parameters. 
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1 Introduction 



The starting of the LHC operation in 2007 brings closer the moment where we will be 
able to see whether the simplest Standard Model (SM) higgs particle really exists PQ, and 
whether the much searched for scale of New Physics (NP) lies in the TeV range [2]. 

Among the many such NP forms, the most widely studied is certainly M = 1 Su- 
persymmetry (SUSY) in four dimensions. Hints supporting the minimal supersymmetric 
model (MSSM), with real electroweak scale parameters, are supplied by the unification 
of the gauge couplings at the scale of ~ 2 • 10 16 GeV and the proximity of this scale 
to the one needed for understanding the neutrino masses and leptogenesis [3]. The easy 
accommodation of Dark Matter in R-parity conserving SUSY models, is also encouraging. 

But, if SUSY candidates are discovered at LHC [U El E], the need to check in detail 
their properties in an international linear collider (PLC) will be overwhelming [8j [9] . 

This note concerns the important option of ILC to run as a 77 collider 7LC 77 , using 
back scattering of e ± -beams by laser photons [101 [Til 12] • In many instances, such a 77 
collider is more efficient for studying new physics, than the standard e _ e + -ILC [131 [HI 

usi rrnj. 

A very important set of processes observable in 7LC 77 , consists of the neutral gauge 
boson production 

77 -> 77, jZ, ZZ , (1) 

which are first realized at the 1-loop level, in both SM and MSSM, and could therefore 
be very sensitive to possible new physics. Their analytical 1-loop study was first per- 
formed for SM in [T71 [THJ, [TJJ], and subsequently extended to include all possible MSSM 
contributions, with real /1 and soft SUSY parameters [20 1 12TI [221 l2~3j- 

It has been found there that, at cm. energies above 200GeV, the dominant helicity 
amplitudes for the processes ([TJ are mostly imaginary and satisfy helicity conservation 
(HC), which forces the sum of the helicities of the two incoming photons to be equal to 
the sum of the helicities of the two outgoing gauge bosons [24]. This striking phenomenon 
mainly comes from the H^-loop, and could make these processes a very interesting field of 
study in JLC 77 . Thus, if e.g. SUSY candidates are discovered at LHC, the study of ([TJ) 
at the contemplated 1% level of accuracy, should provide important checks of their nature. 

In studying ([TJ, the independent helicity amplitudes have been expressed in terms 
of the simplest Veltman-Passarino functions Bo, Co, Do [25]. The resulting formulae are 
quite manageable for the 77 and jZ cases. For 77 — ► ZZ though, the analytic expressions 
of [221 [23], particularly for MSSM, are so complicated, that a public numerical code will 
be very useful. 

The purpose of the present paper is to release such code "gamgamZZ", applying to 
SM and MSSM, for any set of real parameters [26J. For both models, the code calculates 
the independent helicity amplitudes for 77 — > ZZ, as well as the observable cross sections 
in LC 77 , for the case where the polarizations of the final ZZ-pa.ii are summed over. The 
definitions of these amplitudes and cross sections are given in Section 2. In Section 3, we 
discuss the code and give some numerical examples. Section 4 contains the Conclusions. 
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2 The 77 —> ZZ helicity amplitudes and cross sections 

We are interested in the process 

7(Pi7^i)7(P2,A 2 ) Z(p 3 , X 3 )Z(p 4 , A 4 ) , (2) 

where the momenta pj and helicities Xj, of the incoming and outgoing particles are indi- 
cated in parentheses. The generic form of the relevant 1-loop Feynman diagrams is indi- 
cated in Fig{Tk,b. The corresponding helicity amplitudes, denoted as3 Fx 1 x 2 x 3 x 4 (Pz, t, u), 
are expressed in terns of the kinematical variables 



s = s 7J = (pi + p 2 f = 1 ' ( 3 ) 

t=(Pi-P3) 2 = m|-|(l-/3 z cos0) , (4) 

c 

u = (px -Pa) 2 = ra| - -(1 + (3 z cos9) , (5) 



where # is the cm scattering angle, and (3z coincides with the Z-velocity in the ZZ cm. 
frame, provided its positive sign is selected from ([3]). The possible values of (Ai,A 2 ) are 
±1; while those of (A 3 , A 4 ) are(±l,0). 

Bose statistics, combined with the Jacob- Wick (JW) phase conventions |27j, demand 



F Xl x 2 x 3 x 4 (Pz,t,u) = Fx 2 x Mi (Pz,t,u){-l) x ^ , (6) 

F Xl x 2 x 3 x 4 Wz,t,u) = F A2AlA3A4 (/? z ,M)(-l) A3 - A4 , (7) 

Fx 1 A2A3A4 (Pz,t,u) = F Xl x 2 x 4 x 3 {f3z,u,t) , (8) 

while CP invariance implies 

F XlXM ((3 z ,t,u) = F_ Xl ,-x 2 ,-x :i ,-x 4 ((3z,t,u)(-l) x ^ . (9) 



These relations allow the calculation of the complete set of the 36 helicity amplitudes, in 
terms of the 10 basic ones 

F +++ _(/3 z ,t,u) , F ++++ (Pz,t,u), 
F+-++(f3z,t,u) , F + „ 00 (fl z ,t,u) , 
F++oo(f3z,t,u) , F +++0 (P Z , t, u), 
F + _ +0 (p z ,t,u) , , 

F ++ _4f3 z ,t,u) , F ++ _ (p z ,t,u) . (10) 

In [221 [23] these amplitudes have been analytically expressed in terms of the -Bo, Co, -Do 
Passarino-Veltman [25J functions. The code presented bellow is based on these expres- 
sions. As far as SM is concerned, these results agree with those of |19j . except for a trivial 

1 Their sign is related to the sign of the iS-matrix through Sx 1 \ 2 \ 3 \ i = 1 + i(27r) 4 <5(j>4 + pz — pi — 

P2)F\ 1 \ 2 \ 3 \ i . 
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misprint in F + _ +0 , pointed out in [23J. 



As expected on the basis of the helicity conservation (HC) theorem, established in 
for any two body process, the dominant helicity amplitudes in either SM or MSSM, 
should obey 

Ai + A 2 = A 3 + A 4 , (11) 

for s and \t\ much larger than all squared masses in the model. In MSSM, HC is an exact 
asymptotic theorem; while in SM, it has only been proven to 1-loop leading logarithmic 
accuracy [24"] . 

This means that, at asymptotic energies and finite angles, all MSSM amplitudes of the 
set ( fTUj) should vanish, except (F ++ ++, F + _ o), which increase logarithmically 

with energy and are mostly imaginary [22J . 

We have checked numerically that the 1-loop results of [221 EH], indeed obey this rule. 
This is particularly striking for -F++oo> where various contributions, often involving ratios 
of gauge, chargino and charged higgs masses, have been seen to cancel each other exactly 
at asymptotic s, \t\. The validity of HC provides a very nice check of our code. 

HC is essentially true also for SM, the only difference being that the amplitudes that 
violate (ITTj) go asymptotically to small generally non-vanishing constants. 

Backscattering polarized laser photons off polarized e T -beams, would realize an 7LC 77 
collider JTU1II2]- Restricting to the case that the final Z polarizations are not studied, the 
observable cross sections at 7LC 77 are given by 



do^^ZZ) f /3 Z \ ^ 2 2 

a cos \ 1287rs / f— ' 

A3A4 

da^irr-^zz) f Z \ ^ 2 |p ,21 ( 1 o\ 

1 a ~~ Too - / A KV+A3A4 - -r + _A 3 A 4 ' t i6 

a cost) \ lzSns I 

v 7 A3A4 

A3A4 

• J 2^ Re ^ F +- x ^ F -+^A > ( 15 ) 

' 7 A 3 A 4 



dcos9 V 1287rs 



^33(77 ~> gg) f j^ W p r F p* 1 / 1fi X 
^ " (j^JZ^^W^J, (16) 



^23(77 ~» gg) 

<icos# 



A3A4 

(^i)EM^ + AaA 4 n-A 3 Aj, 



(17) 



A3A4 

in terms of the helicity amplitudes satisfying flBl fTOj) 
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3 The gamgamZZ code 



The released code gamgamZZ. tar. gz [26], consists of four FORTRAN codes. Two of them, 
sml_l and mssml_l, calculate the 10 basic 77 — > ZZ helicity amplitudes (flOl) . in SM 
and MSSM respectively. The remaining helicity amplitudes may then be obtained using 
(JHHSJ) - The needed higgs width may be calculated from e.g. [28J. 

For compiling these codes, Loop Tools [29J is required. All input parameters are con- 
tained in input files called "susypar.in" in both, the SM and the MSSM case. All dimen- 
sional parameters are in TeV. The results for each 77 cm. energy, are given for N$ + 1 
equal steps of cos#, starting from cos(9 a ) and ending in cos(0&). The integer parameter 
N g , and the angles 9 a and 9b in degrees, are specified in the in-file. 

As an example we give in Table 1 the results of the sml_l output, for 5 of the 
10 helicity amplitudes listed in (flUj) and N$ = 8 . A similar output is created for the 
remaining five basic amplitudes. The cm. energy is taken at ^/s = 0.5TeV. The exact 
form of all outputs is explained in the Readme.dat file accompanying the codes. 



Table 1: Form of the sml_l output for 5 of the 10 helicity amplitudes listed in ffTQ]) . A 
similar form appears in the same output for the description of the remaining five 

amplitudes of ( jTOl) . 





cos 6 


««*++++ 


ImF ++++ 




ImF + _ ++ 


ReF + _ 00 


ImF + _ 00 


ReF + + 00 


imF ++00 


ReF +++0 


ImF +++0 




0.866 


0.144E-1 


0.162 


0.671E-3 


0.179E-2 


-0.328E-3 


0.360E-2 


-0.267E-2 


-0.473E-2 


0.551E-3 


0.183E-2 




0.707 


0.104E-1 


0.137 


0.655E-3 


0.204E-2 


-0.561E-3 


0.633E-2 


-0.246E-2 


-0.426E-2 


0.351E-3 


0.124E-2 




0.500 


0.772E-2 


0.120 


0.472E-3 


0.228E-2 


-0.832E-3 


0.857E-2 


-0.227E-2 


-0.390E-2 


0.164E-3 


0.783E-3 




0.259 


0.631E-2 


0.110 


0.322E-3 


0.244E-2 


-0.104E-2 


1.000E-2 


-0.215E-2 


-0.368E-2 


0.594E-4 


0.381E-3 




0.0 


0.588E-2 


0.107 


0.2679E-3 


0.250E-2 


-0.112E-2 


0.105E-1 


-0.211E-2 


-0.361E-2 


0.0 


0.0 




-0.259 


0.631E-2 


0.110 


0.322E-3 


0.244E-2 


-0.104E-2 


1.000E-2 


-0.215E-2 


-0.368E-2 


-0.594E-4 


-0.381E-3 




-0.500 


0.771E-2 


0.120 


0.472E-3 


0.228E-2 


-0.832E-3 


0.857E-2 


-0.227E-2 


-0.390E-2 


-0.164E-3 


-0.783E-3 




-0.707 


0.104E-1 


0.137 


0.655E-3 


0.204E-2 


-0.561E-3 


0.633E-2 


-0.246E-2 


-0.426E-2 


-0.351E-3 


-0.124E-2 




-0.866 


0.1443E-1 


0.162 


0.671E-3 


0.179E-2 


-0.328E-3 


0.360E-2 


-0.267E-2 


-0.473E-2 


-0.551E-3 


-0.183E-2 



The other two codes, sm2 and mssm2 are compiled similarly and calculate the differ- 
ential cross sections (11211171) in SM and MSSM respectively. They follow a similar format, 
fully explained in the aforementioned Readme.dat file. These results do not of course 
include the JLC 77 luminosity functions. But it should be easy to combine them with a 
code like CompAZ in order to calculate ZZ production in a realistic iXC 77 environment 

The sml_l, sm2 codes give the exact SM 1-loop contribution from the diagrams in 
FigOKb. 

Correspondingly, the MSSM codes mssmLl, mssm2 contain the exact W, higgs, 
quark and lepton loop contributions due to the box and triangular diagrams in Fig{Tk,b. 
The sfermion and chargino 1-loop triangular contributions from the diagram £0 in Figflb 
are also calculated exactly. 

In these codes though, we have neglected the sfermion and chargino mixing effects, 
arising from the non-diagonal Z couplings, in the box diagrams of Figfjji. This is al- 
ways a very good approximation for the (very small) sfermion box contribution. It is 

2 The chargino or sfermion mixing cannot induce a non-diagonal Z-coupling contribution in this case. 
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also true for the chargino box contribution for all helicity amplitudes, except for the am- 
plitudes involving two longitudinal Z bosons. The reason for this is because the high 
energy chargino-box contribution to the production of a longitudinal gauge boson is de- 
termined by the gaugino-higgsino-higgs coupling [30J, for which chargino mixing cannot 
be neglected. The only amplitudes of (fTUl) belonging to this category are -F ++0 o and F + _ o, 
for which the chargino box contribution is approximated by the asymptotic part of (A. 37, 
A. 38) of [22J and (A. 57, A. 58, A. 63, A. 64 ) of [23] ■ This is an adequate approximation, 
since these contributions are rather energy independent and very small. 

If it becomes necessary, the exact chargino 1-loop box contribution could be included 
in our codes, by using the full analytical expressions appearing in (A.36-A.66) of the Ap- 
pendix of [23] . 

It is also important to note that when using LoopTools [2H] for calculating the needed 
Passarino-Veltman functions, problems may appear related to the small masses of the 
quarks and leptons of the first two families. We have checked that if the LoopToolsl.2 
version is used, there is no difficulty. For LoopTools2 and newer versions though, the 
masses of the quarks and charge leptons of the first two families should be increased to 
around 7 or 8 MeV, to overcome this problem. This does not reduce the accuracy of our 
results, and is easily done by changing the input parameters in the in-file. 

As an example, we present in FigsfSH the differential and integrated cross sections <To 
and (T22, defined in ( ff2l [TBI , for SM (with m# = 130GeV) and for the two benchmark 
MSSM models SPSla' [5] and "light higgs" [6]. SPSla' is an mSUGRA model with 
(mo = 70 GeV, m 1/2 = 250 GeV, A = -300 GeV, tan/3 = 10); while in "light higgs", 
the electroweak scale parameters have been selected as Mi = 150GeV, M 2 = 300GeV, 
M 3 = 600GeV, A t = A b = A T = 700GeV, fj, = 700GeV, m A o = 104 GeV, tan/3 = 34.4 
and all sfermion masses are set at 340GeV [6]. In all cases, the top mass was put at 
m t = 173.8GeV. 

The specific cross sections a , <t 2 2 have been chosen because they often give the largest 
effect. They are also the ones that receive contributions from the s-channejf] higgs ex- 
change diagrams of Figffb. Thus, if the SM higgs is above the ZZ-threshold but not too 
heavy, a higgs peak should be visible in the a and cr 2 2-plots [23]. This applies also to 
MSSM, provided that the heavier neutral higgs particle H° is above the ZZ-threshold 
and its couplings are appropriate; see examples in [23] • Among the scenarios considered 
here though, only SPSla', implying m#o = 0.424TeV, satisfies the 2mz constraint. But 
the H° contribution is so weak in this model, that no peak is visible in FigsJH 

In any case, the SUSY modifications to the SM predictions for cr and CT22, are gen- 
erally energy dependent and sometimes considerably larger than 1%. Therefore, SUSY 
effects in an JLC 77 study of 77 — > ZZ, may be observable and useful. 

Before finishing this Section, it is important to remember that nowadays, and surely 
3 In fact the only cross sections that receive such contributions are <7q, 022, 0-23; compare (|12II17|) . 
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much more in the future, computer codes like FORMCALC [32j or GRACE-LOOP [33] 
are being developed which, starting from a given Lagrangian, construct automatically all 
verticesj, draw the necessary diagrams and calculate all amplitudes in an analytic (or 
semi-analytic) way Applying to 77 — > ZZ, this means that all work contained in e.g. 
[22j [23] is done in an automatic way, obviously deriving similar analytic results. These 
analytic results should then be transformed to a FORTRAN code, with no need of extra 
labor. Under the same choice of approximations, this FORTRAN code should be similar 
to the one presented here, admittedly after considerable labor. In any case though, no 
other such 77 — > ZZ code seems to exist at present. 

4 Conclusions 

We have presented here the code gamgamZZ, contained in gamgamZZ.tar.gz [26]. This 
consists of four codes. Codes sml_l and mssmLl calculate the 10 basic helicity ampli- 
tudes (Tl0|) for 77 — > ZZ, in SM and MSSM respectively; while codes sm2 and mssm2 
calculate the cross sections in (ll2H17p . When the Z polarization is not looked at, these 
cross sections constitute all possible observables in an ILC^. 

The accompanying Readme.dat file contains all necessary instructions for compiling 
and running the codes [26]. The only restriction, applying to the MSSM case, is that the 
soft SUSY breaking parameters and fi must be real. But if interest arises, this restriction 
could be easily overcome, by changing certain couplings in the codes. 

Since the codes sml_l and mssmLl determine all possible helicity amplitudes, the 
results should also be useful for calculating any kind of Z-polarization effects observable 
in JLC 77 , in either SM or MSSM. 

Acknowledgments : 
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4 Possibly employing useful approximations, like neglecting chargino mixing in boxes, as we have done 
here. 
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Figure 2: Differential cross section g?o"o(77 — > ZZ)/d cos 9 as a function of cosd, in SM, 
SPSla' [5] and the light higgs [6] benchmark model. 
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Figure 3: Differential cross section ^022(77 — > ZZ)/d cos 9 as a function of cos#, in SM, 
SPSla' [5] and the light higgs [6] benchmark model. 
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Figure 4: Integrated cross sections 00(77 "~ * an d 022(77 "~ > ZZ), as functions of the 
center of mass photon-photon energy \fs = ^/s 77 , in SM, SPSla' [5] and the light higgs 
[6] benchmark model. 
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